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ABSTRACT 

The traditional method of data presentation for heat exchanger 
surfaces does not permit comparison of individual surface types in 
any simple manner, nor is the published experimental data plentiful 
enough to conduct a meaningful study of the effect of varying 
geometries on surface performance. This data is most commonly 
presented in terms of heat transfer coefficients and friction factors 
referenced to the exposed area as a function of Reynolds number 
based on the minimum free flow area. There are six basic geometric 
parameters for a triangular pitch circular finned tube bank: tube 
diameter; tube transverse and diagonal pitch; fin height, thickness, 
and spacing. 

Baskerville (3) developed a method of surface comparison for 
circular finned tube heat exchangers. Biery (4) developed a method 
for prediction of neat transfer for finned tube banks in crossflow. 
Jameson, et al., (5) developed an empirical relation for correlation 
of pressure drop data across triangular pitch circular finned tube 
Danks. 

The correlations for heat transfer and pressure drop for trian- 
gular pitch circular finned tube banks in crossflow are used with 
the comparison method to study the effect of varying geometry on 
surface performance. Optimal fin spacing and height are found for 
specific fin thicknesses, tube diameters, and tube pitch. Performance 


trends for the six possible geometric parameters are established. 
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NOMENCLATURE 
Dimensionless factor; defined by equation (22b). 


Minimum free flow dimension between tubes measured 
diagonally between adjacent rows from base to base as 


if fins were not present; see Figure 3, [FT]. 


Minimum free flow dimension between tubes measured 
transversely to the direction of flow from base to base 


as if fins were not present; see Figure 3, [FT]. 


Heat transfer area of the bare base tubular surfaces; 


: ; 2 
equals length times heated perimeter [FT]. 
Exchanger minimum flow area for gas: (minimum area 
can be transverse to flow or diagonal to flow; see 


Figure 3). faye. 


Total fin area of a tube of length L; defined by 


equation (A14), feree 

Exchanger flow passage frontal area ignoring” any 
x : 2 

enhancing surfaces; see Figure 3, i egee © 


Ratio of fin area to total transfer area; defined by 


equation (A15). 


Total air heat transfer surface area (same as total 
surface area exposed to gas crossflow), including fin 


surface and edges of fins (ere 


Longitudinal center-to-center tube spacing; see Figure 3, 


[FT]. 





Constant used in the computation of friction factor (f) 
for flow over smooth tube banks; defined by equation 
(32). 


Constant required for the computation of radial fin 


efficiency; defined in Figure 5. 
Dimensionless coefficient defined in equation (23). 
Specific heat at constant pressure, [BTU/LBM-°F]. 


Constant used in calculation of friction factor for gas 


flow over smooth tube banks; defined by equation (31). 


Constant required for the computation of radial fin 


efficiency; defined in Figure 5. 
Equivalent diameter defined by equation (20). 


Equivalent diameter of ae finned tube; defined by 


equation (A118). 
Fin diameter, [FT]. 


Equivalent diameter of a triangular pitch finned tube 


bundle; defined by equation (26). 

Hydraulic diameter; defined by equation (1a), [FT]. 
Ratio of hydraulic diameter to the equivalent diameter 
for an equilateral triangular pitch infinite smooth tube 


bank (ETP-I-STB); defined by equation (21). 


Nominal diameter; defined by equation (1b), [FT]. 





EiP-l-Si6 


ft! 


Tube outside diameter, [FT]. 


Factor defined by equation (20). 


Equilateral triangular pitch infinite smooth tube bank. 


Transverse center-to-center tube spacing; see Figure 3, 


[FT]. 

Friction factor predicted by the Jameson (5) correlation 
and based on total area (AL); defined by equation 
UN 

Friction factor based on total area (A_); defined by 
equation (4a). 


Friction factor used by Jameson (5); defined by equation 


(29). 

Number of fins per inch; defined by equation (46). 

Friction factor based on base area (A); defined by 

equation (4b). 

Friction factor based on total area (A_) for a smooth 

surface, defined by equations (31) and (32). 

32.174 [LBM/LBF] fapceerie 

Mass flux based on minimum free flow. area, At: 
: 2 

defined by equation (2a), [LBM/HR-FT“]. 

Dimensionless geometry factor; defined by equation (A7). 

Mass flux based on free flow area, [Ae defined by 


equation (2b), [LBM/HR-FT°] : 





JBSy 


Heat transfer coefficient based on total area, [A 
2 


defined by equation (5a), [BTU/HR-FT*~--F]. 


Heat transfer coefficient based on base area, [A 
defined by equation (5b), fen Hee Toe F]. 

Projected perimeter of a finned tube; defined by 
equation A117, [FT]. 


Colburn j-factor based on total area, [A defined by 


=} 
equation (7a). 


Colburn j-factor based on outside tube diameter [D J; 


defined by equation (A22). 


Colburn j-factor predicted by Biery's [4] method as 
Seutiimned in Chapter 3. 


Colburn j-factor based on base area, [A defined by 


5d 
equation (7b). 


Colburn j-factor based on experimental data from Kays & 


London [2]. 


Colburn j-factor based on equivalent diameter, [do]: 


defined by equation (A20). 


Colburn j-factor based on total area, [A+], for a smooth 


surface, defined by equation (24). 
Thermal conductivity [BTU/HR-FT-F ]. 


Fin height defined by Real eo for rectangular fin sheets, 


see Figure 5, [FT]. 


HD 


mee 


Sal 


Nu 


NTU 


TFLcorr 


pele 


10 


Flow length along the axis of the tubes (i.e. for the 


fluid inside the tubes). 


Component of fin efficiency ( 
HO) 


” Pe defined by equation 


Number of tube rows in a tube bank. 


Parameter used in algorithm of Figure A4 for calculating 


DL when minimum free flow area occurs in the diagonal 


space. 


Parameter used in algorithmof Figure A4 for calculating 


DE when minimum free flow area occurs in the transverse 


space. 


Parameter used in algorithm of Figure A4 for calculating 


nominal diameter [DJ- 

Symbol used in Kays and London [2] for PR. 
Symbol used in Kays and London [2] for RE. 
Stanton Number [h/G_C,]. 

Nusselt number; defined by equation (6a). 
Nusselt number; defined by equation (6b). 
Number of transfer units. 

Corrected fin efficiency per Figure 5. 


Fin temperature effectiveness of a longitudinal fin of 


rectangular profile; defined by (9). 





A 
Teorr 


AP 


Pr 


RE 


a= 


RE 


i 


Correction factor applied to (9) to permit the use 


Tele 
of the simpler [Tanh ml]/ml equation for fin temperature 
effectiveness in lieu of Bessel Function type equations. 


See Figure 4. 
Total surface temperature effectiveness, defined by (8). 
Pumping Power, [HP]. 


Property constant for fin material and air properties; 


defined by equation (A28). 
Seieticnmoresore dnes, (lse7er- 1 
Prandtl Number [HC /K]. 

Heat transfer rate, [BTU/HR]. 


Reynolds number based on minimum free flow area, [A_]; 


€ 
defined by equation (3a) and labeled N, by ref [2]. 


R 
Reynolds number based on the equivalent diameter of a 
finned tube handle, [D J, and the mass flux based on 
the minimum free flow area, [G J; defined by equation 
(28). 


Reynolds number based on free flow area, [A defined 


mJ; 
by equation (3b). 


Radius of enhanced (finned) surface, [FT]. 
Dimensionless geometry factor; defined by equation (A5). 


Hydraulic radius, equals hydraulic diameter divided by 


four for a round tube, [FT]. 





——_ 
a 






i 


Outside radius of base tube, [FT]. 

Ratio of spacing between fins to tube outside diameter. 
Temperature, [°F]. 

Cold gas inlet temperature, TVs Fewle 

Ratio of fin thickness to tube outside diameter, [ 8 /D 1. 
Hot fluid inlet temperature, [°F]. 


Volume, ened 


Volume one side; volume of the gas side of the heat 


exchanger on a per tube basis. (ere 


Total volume of heat exchanger on a= per tube basis 
fee 


Mass flow rate, [LBM/HR]. 
Principle dimensions of heat exchanger, [FT]. 
Ratio of diagonal tube spacing to tube outside diameter. 


Ratio of longitudinal tube spacing to tube outside 


diameter, equals b divided by DS: 
The greater of Xr or Xp- 


Ratio of longitudinal spacing to tube outside diameter 


for an equilateral pitch tube bank. 





: 
“ 








>), 


is 


Ratio of transverse tube spacing to tube outside diameter 


for an equilateral triangular pitch tube bank. 


Miscellaneous 


Subscripts 


= Oe 


MAX 


OS 


In tube geometry, angle defined by the _ longitudinal! 
centerline and the diagonal centerline; see Figure 4, 


[degrees }. 


Ratio of total transfer area on one side of the exchanger 


to total volume of the exchanger, labeled as a in reference 
2. S 
(2a AUR A= ual 


Fin thickness, or average fin thickness if tapered, [FT]. 
Viscosity, [LBM/HR-FT]. 


Density, [LBM/FT?] : 


Diagonal! 

Enhanced surface 

Equilateral triangular pitch. 

The maximum value a variable may assume. 
Heat exchanger fin material 

Outside 

One side 


Predicted 
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Re /> Reference value for indicated parameter. 

S Smoth surface 

iF Total 

1 Tube side; hotter fluid flowing inside smooth tubes. 

Z Finned side, colder fluid flowing over’ finned tube 


banks. 
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eimMtrOGiuiectron 
A. Purpose 

The objective of any heat exchanger comparison is to enable 
the designer to select that unique surface which can be termed as 


the most beneficial to his particular application. The comparison 


technique should be logical, easily implemented, and inherently 
aeecurate. 
Enhanced surfaces, in the form of plate fins and _ finned 


tubes, have long been used in heat exchanger design to improve 
heat transfer characteristics. Many have presented methods. of 
comparing the performance of various types of surface designs. 

Soland, et al., (1) developed a method of comparison that 
permits performance comparisons of all the Kays and London 2) 
plate finned and unfinned surfaces. Baskerville (3) developed a 
method of comparison, using theory and principle inherent in 
Soland's method as a basis, that permits performance comparisons 
of crossflow finned tube surfaces presented in Kays and London. 
Biery (4) developed a method to predict heat transfer coefficients 
for gas flow normal to smooth and finned tubes. Jameson, et al (5) 
developed an empirical relation that correlates pressure drop data 
for circular finned tubes. 

The purpose of this paper is to: 

1. Employ the methods and correlations developed by Biery, 

Jameson, and Baskerville to predict the performance of a 


crossflow circular finned tube surface for any given 


geometry. 








Zh 


2. Evaluate the accuracy of the predicted performance by 
comparison to surfaces sireganies in Kays and London. 

3. Use the performance predictions to evaluate the effect of 
varying geometry on the performance of = crossflow 
circular finned tube surfaces. 

4. Develop graphs and tables that will assist the designer 
in selecting the geometry of a crossflow circular finned 
tube surface which will be the most beneficial for his 
application. 

B. Background 

Crossflow finned tubular heat exchangers have found 
applications as gas turbine plant intercoolers, aircraft engine 
and electronic coolers, air conditioning units, and various other 
heat exchanger uses. These applications most often involve 
gas-liquid service. [See Figure 1]. 

The performance of a given heat exchanger in which the 
ratio of the heat transfer coefficients between the two working 
fluids is greater than about three is markedly improved by the 
employment of enhanced surfaces. The increased surface area 
inherent with the employment of fins yields a net improvement in 
heat transfer despite lower heat transfer coefficients. 

Optimum design is often used synonomously with "lowest 
cost"'. Cost is frequently defined in terms of amortized acquisition 
cost, operating cost, or impact cost in terms of either weight or 
volume. An example of the latter might be aircraft engine 
component design where weight and volume are valuable 


quantities which can be directly equated to dollar cost. 


we a 
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Estimation of dollar cost of a particular heat exchanger 
design is difficult, primarily due to the proprietary nature of the 
required cost estimating relations. The term "cost" or "optimum", 
therefore, as used in this paper, will most generally refer to 
"impact cost'’ and will serve as a means of relative ranking of 
proposed surfaces specifically by required volume, pumping power 


and heat transfer (NTU) comparison. 
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ll. Comparison Method 

Traditionally, data is presented in terms of heat transfer 
Colburn modulus, j, and friction factor, f, based on the total 
exposed area, At; as a function of Reynolds number, RE, based on 
minimum free flow area, Acs and a hydraulic diameter, DL» of the 
flow passage. [See Figure 2] 

Baskerville's (3) comparison method converts these j and f 
magnitudes to new quantities Jn and uN based on the heat transfer 


area of the bare base tube surface, AL? and a Reynolds number, 
RE) which is referenced to the open flow passage area, Ae; as 
though fins were not present. The effect of the fins is accounted for 
as an increased heat flux and hence larger h for the bare base 
tube surface area. In order to incorporate the effect of the fins 
Into Jnp the metal conductivity must be specified. 

To convert data presented in the format of Kays and London 
(2) to the new basis applicable to the crossflow finned tube case, 
Baskerville (3) derived various ratios from basic heat transfer 
definitions and theory presented by Soland (1). Important assump- 
tions are that the heat transfer resistance of the tube walls is 
negligible and that the controlling heat transfer resistance is 
assumed to be on the finned side of the heat exchanger. Comparison 
assumes that the flow rate, W, the hot fluid inlet temperature, 
eter ana © the cold gas inlet temperature, Te ay are held 
constant. Basic definitions are presented in Table |. Ratios for the 
two cases of minimum free flow area which may occur [See Figures 3 


and 4] are presented in Table Il. Fin efficiency is calculated per 


RIGure oO. 
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FIGURE 2: Sample data plot showing traditional method of data 


Presentation. 
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FIG. 92 


PINNED CIRCULAR TURES 
SURFACE CP = 1.34 


Tube cateide diameter - 0.38 in. 

Fin piteh - 7.54 per inch 

Flow paseage hydraulic diameter - 4, 00.0154 ft. 
Pin thickmeee (average)® = 0.018 i: 

Free-flow area/frontel eres - ¢ 00.558 

Heat transfer ares/total volume — x #140 ft, eee3 
Pin area/total area - 0.892 


Notes Experimental uncertaioty for heat transfer 
reeulte poseibly eomewhat greater than the 
nominal +5% quoted for the other surfacee 
because of the neceesity of eatimating a 
contact reeietance in the bi-metal tubee. 


@— Pings slightly tapered. 
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CASE [|: Minimum Free Flow Area Occurs 


In Spaces Transverse To Flow 


*Minimum Free Flow 


Dimension 








lime att 
FLOW 
A, =al C+ L 
b 
ae = 7 DOL oy = W/A. = W/a,L x, Do 
4 7 
ee i Aeb _ 4a,b _ 4D x [X,-1] 
a% N A 7D 7 
b Oo 

CASE [Il: Minimum Free Flow Area Occurs b 


< 
In Diagonal Spaces 4 
*Minimum Free Flow q 
imensi - 
Dimension s 
D 
O 


FLOW 


<a 
II 


FbL 
5 _4Acb . 8ayb g 8D xX, [Xp)-1] = 
NA 7D rg 

b fo) 


FIGURE 3: The Two Cases of Minimum Free Flow Area Occurence and 








the Associated Geometric Relations of Each Case. 
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CASE |: Minimum Free flow area Transverse to Flow Direction 


Total Volume on one Side: 


7 2 
= [XD] LX, DJL =r Oy L/4 


= dimension into paper along flow length of tubes. 


c Volume: VEZ, 


= [XD . = D,] L XP D, 


CASE IIl: Minimum Free Flow Area in the Diagonals 
ae oe ee. ee <— XD, 


Total Volume on one Side: C ) * 


sina- 7D 2 Ds 
re) 
2} 


D 


vARb" Volume: 





FIGURE 4: Calculation of [ALB ]/Vo. for both Case | and Il. 





Compute fin efficiency assuming a 
Longitudinal fin of Rectangular 
Profile: 
Tanh nl 
Veen aT 


590 


where: ; "a 
ne ae 


What geometric shape describas the 
fin being employed on the tube 
surface? 






Radial Fin of Rectangular 
Profile 


Tapered (Hyperbolic) Rad- Rectangular Profile fin 
jal fin sheet 
{see below for | defn] 





4n 4n 0 &n 


Ss a) BS 
correction correction correction 


C, = 3, in(n Ly) 
Rectangular Profile 
whereas: using r (staggered bank) 


50 
rem ln 2 = 


D,= ~.02426 - .4521 la ‘ 
J = ae 


Vtecorrected - 


4 .orrection 











Compute Surface n: 


A 
a _ . ae "t-corrected! 


FIGURE 5: Computing fin efficiency of a Radial fin, Ref (3). 
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Further, Baskerville (3) shows that for any crossflow heat 
exchanger, where for the same fluid at the same _ temperature, 
u and p are constant, the fluid pumping power per unit volume 


on one side is given by: 


3 c 


P W AP 2 eRe 3 
ij 
: _ NUON Azb ; FREY Aeb (17) 
Viet dy 

“5g “58 ‘ So? DN Yo DN Mos 
And for the flow rate, W, also constant: 
NTU ; mela _ Au jyReyAeb ; jyReyAeb (18) 
V V__WC pr2/ sy pe p * “os 
os os Pp n os n 

Using equations (17) and (18), it is possible to construct a 


performance parameteter plot of NTU/V 56 versus P/Voe from which 
useful performance comparisons between two heat exchangers for four 
different criteria are obtained. The four’ different comparisons, 
depicted in Figure 6, predict relative performance in the following 
categories: 

Case A: Same heat exchanger shape and volume. 

Case B: Same volume and pumping power. 

Case C: Same pumping power and number of transfer units. 

Case D: Same volume and number of transfer units. 

Baskerville (3) plotted performance curves of all the finned 
surfaces of Kays and London (2). Figure 7 shows the resulting per- 
formance parameter plot presenting only the highest performance 


curve at each nominal diameter (D,))- The study considered sixteen 
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FIGURE 6: Performance Parameter Curves For Two Hypothetical 
Surfaces Showing Points Used in Sample Comparisons. 
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FIGURE 7 


PERFORMANCE PARAMETER CURVE FOR FINNED TUBULAR SURFACES 


3000 














Si) 


different surfaces having eleven different nominal diameters. Table 
lll lists all of the Kays and London surfaces considered by type 
and surface designation in addition to the one surface supplied by 
the Trane Corporation. 

Baskerville (3) also offered a procedure which can be used 
for sizing crossflow finned tubular heat exchangers. Two designs of 
an intercooler for a 5000 SHP gas turbine plant were undertaken 
using surfaces 98C and 94. Table IV lists the fluid prcperties 
assumed for both designs. The principle dimensions using surface 


98C and then surface 94 were calculated to be: 


Surface 98C: Surface 94: 

X=6.43 ft SS eon et 

Meze.11 ft Y=10.89 ft 

2-063. Tt ZN FOr tt 

Total Volume = 105.46 ft? Total Volume = 28.68 ae 


Both designs satisfy the thermodynamic and geometric con- 
straints. However, use of surface 94 results in a 73% reduction in 
volume, and is seen from Figure 7 to be notably superior to surface 
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TABLE IV 


5000 SHP Gas Turbine Plant Intercooler Fluid Properties 


Properties Side | Side || 
Faken Water Dry Air 
W 400, 000 lbm/hr 200, 667 lbm/hr 
ie: 60° F WE” TE 
in 
i a? ele 80° F 
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A performance parameter plot is constructed (Figure 8) for 
surfaces 98C and 94 to quantitatively predict the volume savings. 
This comparison corresponds to Case C. From the design point of the 
reference heat exhchanger (surface 98C) a line of slope 1 is 
constructed. The intersection of this line with the surface 94 will 
yield the performance parameters at the design point of the new 
heat exchanger. The predicted volume reduction is found from the 
ratio of either the ordinates or abscissas of the two surfaces. 

The performance parameters of the reference heat exchanger 


design (surface 98C) are: (See Figure 8) 








Niue = 401.20 (1/in~) 
V 

OSs 

P = 2.93 10!7 (pare 
v 

OSs 


The performance parameters that would result using surface 94 


are: (See Figure 8) 








Nails 1634.68) (1/in 

V 

OSs 

Be = - 1.19 WoL faa 
V 

Os 


Since pumping power and number of transfer units are the 


same for both surfaces it is seen that: 


Ba p | 
Vi. | 98C Vo. | 98c Meets: 
a, ae - "a a, es | Vos| 

V 94 V..| 94 98C 


SAS. 


OS 


Therefore, a predicted volume reduction of 75% is anticipated 


using the comparison method proposed by Baskerville. This compares 





She, 





P/V. Caine) 


FIGURE 8 


Performance Parameter Plot for the two surfaces used in the 


Design Example. 
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quite favorably with the heat exchanger design calculations made 
for these surfaces and therefore shows that the performance 
parameters for a surface allow comparison of surfaces without going 
through a complete heat exchanger design calculation. Where there 
may well be hundreds of different surfaces to be considered, such 
as the parametric study of geometry effects on crossflow circular 
finned tube heat exchangers conducted in this paper, the value of 


such a comparison technique cannot be overestimated. 
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lll. Heat Transfer Correlations 

Biery (4) presents a method to predict heat transfer 
coefficients for gas flow normal to finned tube banks with 
triangular pitch. A transformation from the actual tube bank to an 
equivalent equilateral triangular pitch infinite smooth tube bank 
(ETP-1-STB) is made and the heat transfer coefficient for this 
ETP-iI-STB is that desired for the real tube bank. 

The correlation for the desired heat transfer Colburn modulus, 


j, its given as: 


273 ee 


j = NooNep = (Guns SY, 


where RE is based upon the hydraulic diameter of Kays and London 
(2) and is given by: 


D. = 4r, = 4A b (la) 





h h e 

oT 
The minimum flow area, Aa often will be in the transverse 
direction to flow. However, for tube banks with aie small x 
compared to Xo the minimum flow area may well be along the 


diagonal of the triangular pitched bank. Expressions for Ac and DE 


are developed in Appendix | for both cases. 
Ci is a function of the transverse distance between tubes, 
XrD 5 and an equivalent diameter, De» Given by 
= rf x 20 
D. DL +E (D. D,) (20a) 


where for S_-D <.04775 inches 
eis Te 


Sy) (20b ) 


and for .04775 inches =S$ -D, = 4 inches 


B 


Bodin OOvi— ree Spe), (20c) 
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and for “peo, 7 4 inches 


E = 0 (20d) 
where 1 is the spacing between fins in inches. Since E is equal 
fo mone Tor Spe 9, less than .04775 inches and equal to zero for 
Spe0, greater than .4 inches, De has the characteristic of being 
equal to the fin diameter, D., for fins with narrow soacing and 
equal to the bare tube diameter, Do for tubes with either no fins 
or with fins with a wide spacing. 


For an ETP-I-STB the equivalent diameter is simply the 


outside diameter of the tube. Therefore, the ratio of its hydraulic 


diameter, D> to the equivalent diameter is: 
(3) _ 4U(X2)q -1] (Xd, ; 4[(X2), -1] ( 4/3/2) (X-), (21) 
» a 7 T 


where (X7) 4 denotes the transverse tube pitch (ratio of transverse 
tube spacing to tube outside diameter) for an ETP-I-STB. To 
transform from the actual tube bank to an equivalent ETP-I-STB, 
the ratio of the hydraulic diameter to the equivalent diameter 


[D, /D. | is equated to the ratio of the hydraulic diameter to the 


equivalent diameter of the ETP-I-STB [(D, /D.), I. (X7), is then 

solved to yield: 

1 
a 2 oe 
(X>), 1 + (4A + 1) (22a) 
Z 
where 

LGN. 
ge Ph (226) 


Or 4 (4/ 3/2) 


Biery found that with this transformation, Ci. was a function 


). and had the property of bringing all of the various 


of only (xX A 


T 
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configurations onto one curve. Correlating the data given by Kays 
and London (2), Jameson (5), Fraas and Ozisik (6), Mirkovic (7), 
and Grimeson (8) with empirical curve fitting resulted in: 


E = .2818(X 


: - .1282 - 0.08263 [1.2 - (X_),] 


| (X) 4 ~ .8924] 


Biery's transformation and correlation for heat transfer 


Ts (23) 


prediction of circular finned tubes is used later in this paper to 
study the geometry effects on the sizing of crossflow circular finned 
tube heat exchangers. It should be noted that the data from which 
this correlation is derived is based on a limited range of tube 
spacing, tube diameters, fin spacing, fin thickness, and fin height. 
Application of this correlation beyond this limited range should be 
recognized as an extrapolation of the empirical data. 

A study was made to compare Biery's method for prediction of 
heat transfer coefficients with other correlations and the data 
presented in Kays and London (2) for circular finned tubes. The 
details of this study are presented in Appendix II. In general it 
was found that Biery's method resulted in heat transfer coefficients 
that agreed with Kays and London data within 20% of the mean 
value with a mean 7% below the Kays and London data and a 
standard deviation of 8%. 

Zhukauskas (9) developed a correlation for heat transfer for 
turbulent flow over staggered tube banks with the appropriate 
nominal diameter, Dy which is also used later in this paper to 
study geometry effects on the sizing of crossflow circular finned 


tube heat exchangers. This correlation for smooth tube bank Colburn 
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Modulus, es is given as: 


.400 
-0.396 ix ieee (D./D,,) 


or (24) 
ee ooo (Noe) 


To justify the use of the above relation it is noted that 
Baskerville (3) used this relation to predict the heat transfer for 
crossflow over bare tube banks presented in Kays and London (2). 
The results showed a maximum deviation between the experimental 
value of up: from Kays and London (2) and the calculated value of 


j, from equation (24) of 7.5%. These results are shown in Figure 9. 
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[V. Pressure Drop Correlations 
In general the friction factor is a function of the Reynolds 


number: 


5 
Cre (25) 


where C, and C., are constants. 

Jameson (5) developed the following empirical relation for an 
"equivalent diameter", D), of a triangular pitch finned tube bundle 
which correlated his data for pressure drop across tube banks in 


which the minimum free flow area is in the space between tubes 


transverse to the air flow. 








D = we (26) 
j 4 a SN 
H pone XD 
a 1 ieee: 5 ee oe _] 
ae e oe 
where qd. is defined by equation (A118). 
Using this expression for dD), equation (25) becomes: 
f=) 15532 RE nae (27) 
J D 
u 
where RED is the Reynolds number based on D, and Ge: 
J 
G_D (28) 
ame 
Rite = 
D, fA 


The friction factor is defined by Jameson (5) for N tube rows 


as follows: 
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This definition differs from the Kays and London (2) definition 


which is given by: 





2 2 
fG N A 
— _ le 2?G N K, 0, (30) 
if 2p c Avie g D 
O C 0 So H 
Equation (30) is the same as equation (4a), except that 
equation (4a) applies to a single row of tubes, where it is 


understood that the pressure drop is directly proportional to the 
number of tube rows. 

To compare the friction factor used in the Jameson correlation, 
Fy, with the friction factor given in the data by Kays and London, 
the geometry factor, (xX, B/D), must be included and the Reynolds 
number converted to a common diameter. This is done in Appendix 
lll. Im general, it was found that the Jameson correlation predicted 
the Kays and ee data to within approximately 20% from the 
mean value with a mean of +7.8% and a standard deviation of 7.3%. 
As noted by Jameson, this correlation does not accurately predict 
the friction factor for finned tubes in which the minimum free flow 
area occurs in the diagonal of the triangular pitch tube bank. Not 
withstanding this difficulty, Gianolio and Cute (11) in their set of 
experiments for pressure-drop measurements across equilateral trian- 
gular pitch finned tubes banks in induced draft, noted the Jameson 
correlation compared the most favorably with their results showing 
a maximum deviation of only 15% compared to their data. In a 
study of the geometry effects on triangular pitch finned tube banks, 
the effect of tube pitch can be demonstrated with smooth tubes. The 


other geometry effects can then accurately be demonstrated by 
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selecting a tube pitch in which the minimum free flow area does not 
occur in the diagonal. Therefore, the Jameson correlation is used 
later to study geometry effects by predicting performance parameters 
for tube banks in which the minimum free flow area occurs in the 
spaces transverse to the air flow. 


To predict the friction factor for smooth tubes, (f Basker- 


s) 
ville (3) derived the following empirical correlations based on the 
Kays and London (2) data for both cases of minimum free flow 
area. 


For Case | where the minimum free flow area occurs in the 


spaces transverse to flow: 
c 


= -0.2407 ~tOs v 
fo = .7184 RE [X.-1] LX, | (31) 
where: 
= = ee 
CL = .828 - 16.596 RE, 
For Case !! where the minimum free flow area occurs in the 
diagonal spaces: 
C 
. -0.1892 ©2703 D 
fg = RE Dxceil) [X, ] (32) 
where: 
= 0.7284 
Cy = 130.1 RE, 
To justify the use of the above relations, it is noted that 


Baskerville used these relations to predict the actual experimental 
friction factor data for flow over bare tube banks as presented in 
Kays and London (2). The results are shown in Figure 10. The 
maximum deviation between the Kays and London (2) experimental 
value of fg and the calculated value from either (31) or (32) is 


approximately 3.0% from the mean value. 
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V. Comparison of Performance Parameter Prediction Program. 
Appendix IV presents an algorithm and calculator program 
that employs the comparison method introduced in Chapter I! and 
the heat transfer and pressure drop correlations introduced in 
Chapters Ill and IV to predict the performance parameters of a 
surface. This program was used to predict the performance 
parameters for the triangular pitch circular finned surfaces in Kays 
and London (2). These performance parameters were then compared 
to the performance parameters calculated by Baskerville (3) from 
the actual heat transfer and pressure drop data given in Kays and 
London (2) for these surfaces. The predicted performance parameters 
and the performance parameters calculated by Baskerville are shown 
graphically in Figure Il, which offers a visual comparison between 
the predicted performance parameters and the performance para- 
meters generated using the Kays and London (2) experimental data. 
Each surface was compared on the basis of same pumping power and 
Number of transfer units (Case C). With this comparison the ratio 
(V__) AV__)._ can be found. The percent deviation in (V__)_ from 


OS Pp OS ¢ Os ££ 


Ve can then be found from this ratio and calculated according 


Tomo qucitionn (33) 


(oe = (33) 
Percent Deviation = cai x 100% 
os e 


The average deviation of the predicted volume one side Wael 


from the calculated (from experimental data) volume one side Was 
is given by Table V. if surfaces 93, 94, and 98D are discounted 


due to the difficulties already cited for these surfaces, the average 


deviation for the remaining surfaces is +11.7% with a standard 
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Figure 11. Comparison of predicted performance parameters’ with 


performance parameters calculated by Baskerville (3) 


from Kays and London (2) experimental data. 
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TABEER SY 
Deviation of performance parameter program of Figure A5 for volume 
one side in predicting experimentally obtained value as presented 


in Kays and London (2) and calculated by Baskerville (3). 


Surface Average Deviation 
92 +23% 
93 +86% 
94 +519, 
96 +249, 
97A +30% 
97B +20% 
98A +5% 
98B -6% 
98C Pie. 
98D +25% 
98E -2% 
99A +24% 


99B +10% 
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deviation of 13.8%. This degree of accuracy should provide the 


correct indications in the parametric study of geometry effects on 


triangular pitch circular’ finned tubes presented in the next 


chapter. 
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ee cometny Eriects on Irlanguilar Pitch Circular uminnedgiwibe Heat 

Exchangers in Crossflow. 
A. General 

A parametric study was conducted using the program of Figure 
AS to determine the effect of varying tube pitch, diameter, and fin 
geometry on exchanger effectiveness. Copper fin material and dry 
air at 90° F on the enhanced surface side was assumed. Each 
surface was compared to a reference surface for the same pumping 
power and number of transfer units. This corresponds to Case C 
discussed in Chapter II. Therefore, the relative figure of merit for 
the surface in question is the ratio of the volume one side of the 
reference surface to the volume one side of the surface in question, 
(V og REE’ Mo: With the volume of the reference surface in the 
numerator, this relative figure of merit can be thought of as the 
number of times larger (or smaller) the heat exchanger composed of 
the reference surface would have to be in order to have the same 
performance as the surface to which it is being compared. The 
reference surface is taken to be a bare tube surface with a one 
inch tube diameter and an equilateral triangular tube pitch of two. 
The reference surface is completely an arbitrary choice. 
B. Effects of Tube Diameter and Tube Pitch 

The performance parameters produced by the program of Figure 
A5 for smooth tubes (no enhancements) were used to study the effect 
of varying tube diameter and tube pitch. Table VI lists the values 


ey which were generated for various values of tube 


os REF! “os 


diameter and tube pitch. These values are plotted on a log-log 


scale in Figure 12. 
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RaBee Vi; (Vos REF/ Ye for various values of tube diameter and 
tube pitch for smooth tubes. The reference surface is a 
smooth tube with a one inch tube diameter and equila- 
teral triangular pitch of two. 

A. (V pee/Vo~ for varied tube diameter and several tube pitches. 

X7= 1, K-=2 AS X7=3.9, 
D (Inches) Xerp=2 Xerp34 Kale Xe X)=2.0 
4] IBS| fe 4.65 316 Sie SieOs 
D274 SVS - = - 
2.74 440 32.4 Taos 2.99* 
Ride) ey BARS = = = 
1.0 16.8, 5 SNS, T2202 Zao" 1.09* 
2.0 7306 10575 - 2 = 
2.5 = = - Bios Palate 
4.0 134 0208 - = = 
BO 2 = = ioe elOar 
1070 .0354 .00545 - a z= 


8. Variation 


D T 


in equilateral triangular tube pitch. 





os b=" D = +49 
1) ees 106 
1.3 4,83 S20 
1.5 2.61 19.3 
2.0 1.0 7.48 
= 0 REA, 249 
Bn, .0918 705 
10.0 sO 135 
* Indicates minimum free flow area occurs in the diagonal! 
(ye X-_+1). 
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CASE 1 ——— 


CASE 2 —--—- 


versus Xr for a bare tube surface for two 


values of diagonal tube piich. 
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Table VI and Figure 12 show the tube diameter and the tube 
pitch for Case 1, minimum free flow area in the transverse space, 


to be approximately proportional! to exchanger volume as follows: 


errs (34) 
Wasi OG 
eee 22 | Bs ps), Coe) 5710 
OS iT : or are eer 
Zs 
Vosg @ Xp > 1.55X, 53.0 (36) 
Vea a Xx cad lox = 370 (37) 
OS Sips : ear A Cr 
This leads to the following relation for Case !, minimum free 
flow area in the transverse space (X_- 2X 5-1) 
A135 Z2\5 inas 
V X 
OS = a5 Xb x Te soe 0, (38) 
Vos) REF (XT) pep XD) REF OS REF 
Case I!1, minimum free flow area in the diagonal (2D, = X>+1), 


is not as well defined. Again Table VI and Figure 12 show for Case 
1 | that the tube diameter and the diagonal tube pitch is 


approximately proportional to exchanger volume as follows: 


Vos * Dy (39) 
3.85 

Vow © Xp 1.1%, <1.5 (40) 
2.5 Ges 

Ue * Xp 1.5 Xp 3.0 (41) 


However, Figure 12 shows that exchanger volume even on a 
log-log scale is a non-linear function of the transverse pitch. 
Therefore, variation in volume versus transverse tube pitch on a 
log-log scale does not lend to a linear approximation even for a 
very limited range. However, since the change in volume is small 


over the possible range of transverse pitch to a first approximation 
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volume can be estimated to be independent of transverse tube pitch 


for Case Il. This leads to the following approximations for Case ||: 
as 
V X 85 D |. (42) 
os es D O ; amt hi 
VosRer |p) rer (OO REF 
225 ieee 
V x : 
22 ~ D 26 | 1:5 = eeieRO 2) 
(Vos) REF (Xp eee (Do REF 
It should be recalled that Case I! occurs when 2X. = Xr+1, 


but that 2X = X7 is a geometric impossibility. As the latter limit is 
approached the tube bank approaches an in-line tube bank. Since 


x must be at least 4} tube diameter, this further restricts the 


range for Case ||. Since: 
2 oe s Z 
xX /4 + XP = Xo (44) 
Case Il is restricted to the following range: 
2 z 
(X+ ap) SER Gal (45) 


This range of possible values is less than a half tube diameter. 
Due to this limited range of geometries in which Case I|1 may occur, 
the minimum free flow area usually occurs in the transverse space 
for most triangular pitch tube banks. 

Finally, a note should be made concerning the discontinuity 
apparent in Figures 12C and 12D for the transition between Case | 
and Case II. The curves for each Case are produced by a separate 
set of correlations for heat transfer and pressure drop. Therefore, 
unless both correlations are exactly correct, a perfect match would 
not be expected. The amount of mismatch is in a sense a measure 
of the total error of the four correlations taken as a group at the 


point of transition between Case | and Case II. 
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C. Effects of Fin Geometry 


The performance parameters produced Dye the prodmam. or efi aqure 
AS for an enhanced surface with circular fins and equilateral 
triangular pitch were used to study the effects of varying fin 
geometry. Table VII lists the values of (Vg) Ree/ Ve generated for 
various values of fin thickness, fin height, and fin spacing for two 
tube diameters (.4 inch and 1.0 inch) and two values of equilateral 
Eelangular pitch (Xetp = 2 and Xertp = 3 Je ihesce Values sare ploued 
in Figure 13. Again, (V og REE is taken to be a smooth tube with a 
one inch diameter and an equilateral triangular pitch of two. Fin 
spacing is given in terms of fin thickness (Sue/te). The greater the 
value of (Vg REE Mos? the smaller is the V _ of the surface in 
question for a given Eee of transfer units and pumping power. 

Figures 13A, 13C, 13E, and 13G show the effect of varying fin 
spacing with all other parameters held constant for several values 


of fin thicknesses. The fin height is taken to be the maximum 


= 4(X ~ 1). An optional value of fin 


possible value, (H_) ETP 


MAX 
spacing, i.e., a maximum value of (V epee Vos and, therefore, a 
minimum value of ee for the surface in question, occurs for each 
case. This optional value of fin spacing ranges from .8 to 1.2 fin 
thicknesses. The higher values of optional fin spacing in terms of 
fin thicknesses corresponds to the thinner fins. 

Figures 13A, 13C, 13E, and 13G show for a given fin spacing 
in terms of fin thicknesses the surface performance is Detter for 


thin fins than for thick fins. This is not necessarily true if fin 


spacing is not expressed in terms of fin thicknesses. Figure 13l 
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plots (V )pee/Vo. versus fin thickness (t 


OE ) for two constant values 


f 
of fin spacing (5S ; 

BF). Again optimal values are found but here 
thicker fins are seen to produce a surface with better performance 
if the fin thickness is less than the optimal fin thickness. Figure 
ieuy plots (Vos REE! Vos versus aan for two values of fin thickness. 
Again it can be seen that for low values of iS oe the thicker fin 
produces the better performing surface. 

The information plotted in Figures 13A, I3@ A VGEs ands) 3Grcean 
be displayed in the form of Figure 13J, since the number of fins 
per inch is given by: 

1 1 (46) 
Dap oe Shape! Goalie}. 

There may also be other possible forms of displaying the same 
information. 

Figures 13B, 13D, 13F, and 13H show the effect of varying fin 
height from the maximum possible value for values of fin spacing 
which were found to be optimal in Figures 13A, 13C, 13E, and 13G, 
respectively. An optimal value of fin height occurs ranging from 
40% to 90% of the maximum possible fin height. For thinner fins the 
opti mal value of fin height occurs at lower fin heights and is more 
pronounced for thinner fins. 

For some optimal values displayed in Figure 13, the optimal 
regions are quite flat. This suggests that in the actual design of a 


surface great savings of fin material might be made by only 


slightly compromising surface performances. 
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TABLE VII. (Vos REE Mos for various values of fin thickness, fin 
height, and fin spacing for two tube diameters (.4 inch 
and 1.0 inch) and two values of equilateral triangular 
pitch (Xetp = 2 and X =a Peas 


| Ete ye! EE 
tube with a one inch tube diameter and equilateral! 


is a bare 


triangular pitch of two. 


A. (Vos REF/ Moe versus fin spacing for four fin thicknesses 


(t, = .005, .01, .02, .04) for an enhanced tube surface 
with D = 1", Xeqp = 2, and H. = (He) ax = 292 
(Vos REF! os jolie 
Spe/te t,=.005 t.=.01 t.=.02 t.=.04 
3 22.4 (75s 3 az 9.60 
4 24.6 eee 4.4 (ORS ES: 
.6 26.3 WS) 8: | gle NGeeZ 
Me! 26.8 INS) 4 (ere loess 
8 Di iO 206 (emete 10.80 
0 DD Rei 20 ay (pee 1Oe7Z 
8 743 UGS, 205 [ee lOn38 
1.5 26.9 20.0 14.6 10.07 
Za0 26.0 ee Sa Bh 42: 
252 24.2 0 12.6 8.67 
3.0 22.6 eae (RLbsKs: 8.01 
Se0 ZA les 2 7.60 


*Indicates Maximum Value. 
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TABLE VII. (Cont) 


3h. (Vos REE Vos versus fin height for four fin thicknesses 


ee -005, .01, .02, .04} for an enhanced tube surface with 


oe = 1" and Xerp = 2 and with fin spacing at the opti mal 
value for H. = (He) ax = 5 
SU cetisveret Us io 
| eR ee 
ee i i ie io 
750 Ql et ZNO Fi or 2 10.84 
~45 Zou l Zao 15.6 eae 
.40 30.4 22.4 1Oe2 lite 
WOOD rea, 22295 oo 1076 
. 30 S27 Paps) Jerais) - 
ne 32.95 - 14.8 S203 
. 20 3236 19-3 - - 
alto Pas) ee) - 1025 - 
8 Zoe Zs - A23 
05 13.9 - = = 


*Indicates maximum value. 
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Weeee Vii. (Cont) 


on (Vg REE Moe versus fin spacing for three fin thicknesses 


(te = OlmeO2 a O4). [Or anmenmanced. tbe curiacem itm De 
= (4 = Z = 
» Xertp 2, and H, (He) ray SOK 
(Vos REF! “os uel 
Spete te = Oe) te = 0.2 te = 0.4 
re S529 2730 20.1 
41.0 S056 Pipa 
.95 - _ 22 
jee0) Ley, 30.9% - 
eto 41.8 - = 
| A a8 30.9 SLPS 
ies 6 SS | S(O) Zee 
20 40.0 29.0 20.2 
S.0 56.0 Zon 17.6 
eS: S220 D2 Lowe 


* Indicates maximum value. 
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TABLE VII. (Cont) 


Drs (Vos REE’ Moe versus fin height for three fin thicknesses 


(te = “Ole 2.02, (04) an emhancea tube —sUrriace switm Dy = 

ae ovand Xetp = 2 and with fin spacing at the optimal 
= — 

value for H. (He) pay alee 


Spetp=1.15 Sa-/t.=1.05 Spe/t.=-95 
H. te O1 te=.02 te O4 
50 41.8 30.9 D7 
45 45.1 33.1 DVT 
40 48.4 255° Dea? 
35 51.4 36.6 2 ea 
30 53.8 a7 24.8 
25 Sore 36.9 23.6 
20 54.3 34.7 21.2 
15 49.8 30.1 eae 
=10 eon7 22 ah 1oge 


*=Indicates maximum value. 





is 


feet Vil. (Cont) 


E. (V Yeee/ Vo. versus fin spacing for four fin thicknesses 


os 

(t. = .005, .01, .02, .04) for an enhanced tube surface 
ith = "t = = = 

wi DS ae XeTp 3 and He (He) ay sO 
Velen! Ve ae 


Spp/te t. = .005 i, = 10! t, = .02 ey 0 
40 14.75 ‘lege 9.11 6.92 
60 15.66 12.51 9.74 ase 
ATI = a a 7.38% 
80 15.88 JD 9.90 2 
825 15.88* Ds 9.90* = 
85 15.88 12.75 = = 

1.0 Suz 12.70 9.84 7.26 
loz 15.48 12.50 9.67 7.98 
1.6 dh A 11.92 9.17 6.77 
2e0 13.83 11.26 8.63 6.44 
2.5 12.80 10.46 7.99 6.02 
3.5 11.03 9.06 7.06 5.30 
4.5 9.63 7.96 6.30 4.75 


ve se e 
*tndicates maximum values. 
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TABLE VI. (Cont) 


F. (V os REE! Mog versus fin height for four fin thicknesses 


(t. = .005, .01, .02, .04) for an enhanced tube surface 

with De =e. XeTp = 3 and with fin spacing at the 
ti | | = = 

optimal value for H. (Hed aray 1L2Or 


(Vos REF “os for 


Spp/tp=-825 Sap /t.=.825 Spp/t,=.825 Sgp/tp=-775 
H, t,=.005 t.=.01 tp=.02 t= .04 
ae lerorehs: [PAR tes; 9.90 7.38 
35) 17.46 ~ = = 
00 - = - fissile 
she, | Syg1ehe) 14.99 eA, 1293 
TT es. - - - Tos 
See, - - - 7.94 
170 20.77 ~ - Tan 
265 ~ - ele - 
HOZS ~ - 11.67* - 
. 60 22 39 16.74 11.64 ol 
ee), - 16.96 11.48 - 
90 RIT fe. Oss - = 
145 24.21 eaters, - = 
42 24.36 = = = 
-40 EE oe 16.48 10.04 5270 
er) 24.24 ~ - = 
7o0 Pe oe = = = 
», 7A9) 1355/7 S209 550 Zico 


*Indicates maximum value. 
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TABLE VII. (Cont) 


G. (Vos REE Voc versus fin spacing for three fin thicknesses 





(t. = .01, .02, .04) for an enhanced tube surface with 
en =o XetTp = 3, vand (Hedray = 1203 
Soglineie’ “oe yee 
) Sae/te t-=0.1 t.=.02 t.=.04 
4 PAO TeV [e2) Preis ites 97 
748) OVE: 7297 eS 799 
-92 - ~ ioe 
ie 2351296 eye 13.46 
I) ON - ieee ~ 
eal 24-025 = = 
eZ 23 eo5 18.14 13235 
Pe 25.01 7 a4 i060 
220 Dll tee Or 12.26 
S20 20.71 [pores ) Ols We 
Ce ees 13.14 9.43 


*Iindicates maximum value. 





PaBpiee Vil. (Cont) 
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tale Ue ee ee versus fin height for three fin thicknesses 


(te = .01, .02, .04) for an enhanced tube surface with 
Dd. ae ye a Xetp = 3, and with fin spacing at the optimal 
value for He = (He) pay PE so 
Bee) ee ioe 
Sye/te=1.1 Spe/t,=1.07 Spe/t,=-95 
He te=.01 t.=.02 tp=.04 
i® AAR OV 18.18 [ont 
Zo 50 roe 77 a o0 
ae 28.74 PAE) one 7 
ol SH Ass: CORN (eros 
As. 33.66 PISS TET eins 
ms 35.06 DiS 13). 
4 Senor: 5) Mes 14.08 | 
20 34.87 20,3 eases 
2 Zoe 15.04 Ses 


*Indicates maximum value. 
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TABLE VII. (Cont) 


| « (Vos REE Moe versus fin thickness for two values of fin 
spacing (Spe = .04, .10) for an enhanced tube surface with 


Boe 1M Mi = 2, erie) Fl, SEI = 2B 


os REF! “os er 


te Spr = .04 Sar =a, 10 
-005 14.13 asks 
01 i, ihe Sie 
310 13.42 Scie 
04 [OV R tees S07 
06 S00) Watsis; 
05 7220 Ue) 
10 6.56 6165 
20 4.13 4.88 
40 Tee 8559 





78 


feeee Vil. (Cont) 


ie (Vos Ree Moe versus fins per inch for two values of fin 


thicknesses (t = .01, .02) for an enhanced tube surface 
with D_ = 1" = = = 
» Xerp = 2, and H, (Hed aay are 
(Vos REF Vos VOR 


ae te=.01 te= O02 
2.4 4.25 4.98 
eye, — Siete] 
eye: at — 
Ses) - SIRES, 

(erste: OES - 

1353 - 12.30 

G77 13.24 Sass, 

20.0 14.79 14.62 

25.0 Neig (/s. Bee Ad 

Zo ie9S jHaweech ee 

83.3 io co7 14.62 

40.0 133.95 202 

45.5 207.65 

47.6 20.90* = 

50.0 ZO56> ~ 

Son” 1923 - 

80.0 15.65 - 


*Indicates maximum value. 
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Figure 13A. (V oper os versus fin spacing for four fin thicknesses 
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(Vg) /pEE/ Vos versus fim spacing for three fin thick- 
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Effect of tube diameter and pitch for enhanced tube surfaces 
can also be seen from Figure 13 from the two values of tube 
diameter and the two values of tube pitch for which (Vo pee/V oe 
was plotted for the various fin geometries. The effect of tube 
diameter and tube pitch on exchanger volume for an enhanced 
surface is less than a bare tube surface. If the optimal points in 
Figure 13 are compared for the two tube diameters and two tube 
pitches, the vee iS approximately doubled for a 2.5 increase in 
tube diameter and increased 1.5 times for a 1.5 increase in 
equilateral triangular tube pitch. For a bare tube surface larger 
changes in Vine would be predicted for the same change in tube 
diameter and equilateral triangular tube pitch. 

D. Summary, Conclusions, and Recommendations for Future Studies. 

Correlations for heat transfer and pressure drop for trian- 
gular pitch circular finned tube banks in crossflow developed by 
Biery (4) and Jameson (5) are used with a method for comparison of 
finned tube surfaces developed by Baskerville (3) to predict the 
performance of a surface with a given geometry and to study the 
effect of varying geometry on the performance of enhanced surfaces. 
There are six basic geometric parameters for a triangular pitch 
circular finned tube bank: tube diameter; tube transverse and 
diagonal pitch; fin height, thickness, and spacing. 

lt was found that for the same heat transfer (NTU) and 
pumping power, the volume of a heat exchanger decreases for 
smaller tube diameter and tube pitch. Optimum values of fin 


spacing and fin height were found for _ several specific fin 
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thicknesses, tube diameters, and pitch. Thinner fins were found to 
be more effective in enhancing a surface than thicker fins for the 
same fin spacing in terms of fin thicknesses. 

This study extrapolated correlations based on experimental 
data for heat transfer and pressure drop beyond the range for 
which data exists. Experimental or theoretical work remains to be 
done to either prove the validity or correct any correlations that 
might be used in a similar parametric comparison study for the 
range of the geometric parameters where optimum performance 
appears to occur. 

Future studies on the geometry effects of heat exchangers in 
crossflow might be extended to geometries not limited to triangular 
pitch tubes, circular fins, or circular tubes. With the proper corre- 
lations parametric studies, such as this one, might be performed for 
all the other geometries a crossflow heat exchanger might assume. 
Use of a high-speed computer would allow plotting families of 
performance curves versus the various geometric parameters. 

Finally, performance might be defined in terms other than 
minimum volume for the same number of transfer units (NTU) and 
pumping power. A promising candidate might be minimum exchanger 
weight for the same NTU and pumping power. This would involve 
minimization of the volume of the material (tubes and fins) of an 
exchanger instead of the actual volume occupied by the exchanger 
for a given performance. Another dimension to this problem would 
be the materials from which the exchanger is constructed from. For 


example, although aluminum is lighter, copper has a higher heat 
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conductivity. However, materials are often dictated by = service 
environment. Therefore, optimum performance for minimum volume of 
exchanger material would be a most useful measure to minimize the 
weight of an exchanger. This would insure the most efficient use of 
heat exchanger material on a per weight basis for a_e given 


performance. 
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APPENDIX | 


The purpose of this appendix is to present the derivation of 
several geometric relationships. 
A. Derivation of total heat transfer area on one side, Ar: Tels 

circular finned tubes with fins of rectangular profile. 

Figure Al shows the layout of a circular finned tube with fins 
of rectangular profile. 

The total heat transfer area including fin surface and edges 
of fins is: 


_ 2 Z 
a L [D+ [2(7D, /4 A) /4) + ™D-t.D. -7D.t,0,1/[(Sgett,)D9]] (An) 


Since, 
D. = D, + 2H,D. (A2) 
and, 
DM) _ Dee 2 
De = (D,-D_)(D,+D.) = 4H.°D.” + 4H,D, (A3) 
A, = LD, ; - [2H.(H, +t. + 1)/(Spe + te] (A4) 
Betting: 
Ae H./(Sae + t.) AS) 
Ar = *LD [1 + 2R, (H, een 1) ] (A6) 
and letting: 
G. =e Lil 2R. (H. oo 1) ] (A7) 
A, = LD G, (A8) 


B. Derivation of minimum flow area for gas, A. 
Figure A2 shows tthe triangular layout of a tube bank in 


crossflow. 
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FIGURE Al 


Circular finned tube with fins of rectangular profile. 


t¢Do De 


Spr Do 


FIGURE A2 


Triangular layout of a tube bank. 
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Minimum flow in the transverse space occurs when ae 2a 4; 


while minimum flow occurs in the diagonal space when 2a, a,-. 


Therefore, minimum flow in the transverse space implies XS 2x ~1, 


D 


and minimum flow in the diagonal space implies X. = (X~+1)/2. In 


D V 


either case, Figure A3 shows fins reduce the minimum flow area. 
Therefore, in the calculation of minimum flow area, the fin 


cross sectional area must be subtracted. 


For Case |: Minimum free flow in the transverse space: 
XS 2X,-1 
and, 
a L[X.D. Pee 2t.D HD /(Syett.)D. | 
= LD (X7 - 2R ete - 1) (A9) 
For Case !!: Minimum free flow in the diagonal space: 
Xp = (X41) /2 
and, 
a L[2X,D = 20) = 4t-D HD /(Sp-tt.)D J 
= 2LD 1X, = AR Ma il | (A10) 
C. Derivation of hydraulic diameter, D_- 
For Case |: Minimum free flow in the transverse space: 
<= = 
K+ Xo! 
and, from equations (la), (A6), and (A9) 
D, ; 4(X_-2Rt,-1) x, Do ; 4X) Do (X7-2Ret.-1) (A111) 
a Paarl) = 8 a, 
T 1+2R, Hett + 7G, 
For Case !!1: Minimum free flow in the diagonal space. 


Xo = (X_+1)/2 





sil 


PIGURE As 


Layout of a tube bank with gas flow normal to paper. 
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and, from equations (la), (A6), and (A10) 
D, 8(X,-2Rt.-1)X, Dd : BX, D(X,-2R.t.-1) (Ace 
7[1+2R.(H,+t,+1) ] 7™Ge 


D. Derivation of ratio of total heat transfer area on one side of the 
exchanger to total volume of the exchanger, 6 . 
From equation (A6) and Figure 3: 


A w{1+2R,(H_+t.+1) | 7G (A13) 


B eee = —_——__ 
Vr XTX, Do X7X, Dd. 


E. Derivation of ratio of fin area to total transfer area, A./Ay: 


It 1s seen that: 


L[2l 7(D,°/4) - 7(D.*/4)] + 7D, ] 


Be 
Spe tte 0, 


= 27LD OR (H+ t. + 1 + t ./2H,) OTA) 


Ag 2R -(H. vie 1 + t,/2H,) 


a mniieeecRe (Ht +1) = 
AD Ie? 2a (HH 


ee ls 2 t-/2H_)/G, (A15) 


if 
ho 
A 
ag 
ore 


F. Derivation of the angle define by the longitudinal centerline and 


the diagonal centerline, a ; see Figure 4. 
X_D/2 (A16) 
2 vaGeram ae Bearetanm in / 2X0) 
a xD = r Le 


mee) 
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G. Derivation of the equivalent diameter used by Jameson (5) and 
Mirkovic (7). 


The projected perimeter of a finned tube, Hes is 27 times the 


distance traveled axially along the tube and over the surface of 
the fins from one end of the tube to the other. (See Figure Al.) 
Therefore, 


27bL(2H.D + t,0 +S D ) (A17) 


io f BF o 
| aL (2R. + 1) 
kK ted. + “pe? , f 


The equivalent diameter of a_ finned tube, da used by 
Jameson (5) and Mirkovic (7) is defined as twice the total air heat 
transfer surface area divided by 7 times the projected perimeter of 
a tube. Therefore, 


2A D [1+2R,(H.+t.+1) ] DG. (A18) 


e eu, 2R +1 2R +1 
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APPENDIX 11 

The purpose of this appendix is to compare the method 
developed by Biery (4) as presented in Chapter II! for predicting 
the Colburn Jj-modulus for gas flow normal to triangular pitch 
finned tube banks with the data given by Kays and London (2) and 
with several other correlations. 

Table Al presented the percent deviation of the Colburn j- 
modulus predicted by Biery's method from the experimental value 
given by Kays and London (2) for the corresponding surface over 
the given range of Reynolds numbers, where 


JB OT Skee x 100% (A19) 


J KEL 


Percent Deviation = 


Table Al shows the Colburn j-modulus predicted by Biery's 
method deviates from the experimental values by approximately 40 
percent for surfaces 93 and 94. However, Kays and London notes 
experimental uncertainty for heat transfer for surface 93 due to the 
necessity of estimating a contact resistance in the bi-metal tubes 
which were used. Another source of error for surfaces 92, 93 and 94 
is, according to the values given, Kays and London calculated the 
hydraulic diameter without accounting for the fins as outlined in 
Appendix 1, Section C. However, the values given for hydraulic 
diameter for the other surfaces, which are data taken from Jameson 
(5), suggest the fins were accounted for in calculating the 
hydraulic diameter. Therefore, due to these errors, surfaces 93 and 


94 are discounted. The remaining body of data shows that Biery's 





i f 
Deviation o 


TABLE Al 


the Biery (4) 
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GolBUurn 
d for 
metho 


j-modulus 


in 


a d 
Si Nn Nn resente Kays Nn 
iP nm m d by 

t 

1 | obtai ed V 

y 

i enta 

i experi 

dicting 

Dine 


London (2). 
Surface* RE=1000 
-27 .1% 
. -41.7% 
: -37 . 3%, 
5 -12.0% 
; -16.3% 
7 -11.5% 
978 on 
" +17 .8% 
. +27 .4% 
+6.3% 
se * 12207 
"i -10.1% 
99A ae 
998 


 & 
a 


figure number. 
J 


+6 


+13 


-4 


steal 


-13. 


ihe 


RE=4000 
5% 
.7% 
1% 
. 8% 
8% 
5% 
.6% 
0% 
.0% 
0% 
4% 


4%, 
3% 


tal 


+7 


=i 


-2 


-14 


-5 


RE=7000 
.7% 
3% 
.0% 
4% 
.6% 
1% 
4% 
6% 
.7% 


9% 


576 
.7% 
5% 


-12 


-6 


RE=10000 
2% 
0% 
9% 
1%, 
7% 
0% 
7% 
1% 
my, 
Sie 
0% 
6% 
2% 


don 
Kays and Lon 
the 

to 


(2) 
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method of predicting the Colburn j-modulus agrees with the 
experimental values given by Kays and London (2) to within an 
average per surface of 20% from the mean value with a mean 7% 
below the Kays and London data and with a standard deviation of 
8% . 

The method given by Biery (4) for predicting Colburn 
j-modulus was also compared to correlations given by Briggs and 
Young (10) and Mirkovic (7). Caution should be exercised when 
comparing these correlations as the Reynolds number for each author 
is based on a different diameter. To compare these correlations a 
factor must be included to convert the correlations to a common 
diameter. The correlation given by Briggs and Young is based on 
the outside tube diameter while the correlation given by Mirkovic is 
based on an equivalent diameter, da as presented in equation 
(A18). 

The Reynolds number used by Biery (4) is based on the 
hydraulic diameter given by Kays and London (2) and equation 
(la). 

The correlation given by Mirkovic (7) when converted to a 
hydraulic diameter basis may be expressed as follows: 

2/3 


> Nes Neya 


= .224)1.1547X Dd, ~ Dd, bh ee) eB S ; D 


ie 
Os =o He e (A220) 


Ju 


lt should be noted that 
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BE _ 1 - F. td. (A21) 
MU TO : 

which ts the way Mirkovic originally included this factor in his 
correlation. It appears Mirkovic intended to develop a correlation 
which included fin thickness as a variable in this factor, although 
it is not. However, the Mirkovic correlation is slight!y dependent 
upon fin thickness when expressed in terms of a hydraulic diameter 
due to the (D, /d,) term. 

Biery (4) noted the Mirkovic correlation gave very fine results 
when used to analyze the Mirkovic (7), Kays and London (2), and 
Jameson (5) data. However, this should be expected since the data 
produced the parameters and coefficients for the correlation. Biery 
reported the overall result in comparing the Mirkovic correlation 
with these data gave an average error of +.3 percent with a 
standard deviation of 3.6 percent. However, since the Mirkovic 
correlation is not directly a function of fin thickness, its more 
general validity outside the narrow range of data from which it 
was drawn is to be suspected. Also, the Mirkovic correlation 
includes a term that subtracts the diameter of the tube from the 
longitudinal distance between tube centers. In a triangular pitch 
tube bank, where the transverse tube pitch is large and the 
longitudinal tube pitch is small, the longitudinal distance between 
the tube centers can approach and sometimes be smaller than the 
diameter of the tube. 


The correlation given by Briggs and Young when converted to 


a hydraulic diameter basis may be expressed as follows: 
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As S Rasy S oe. 0 319 Rew 39 (A22) 
i = NIN | ae a eee 
JBEY ST Pr t H D 
f i O 
Gianolio and Cuti (11) confirmed in a set of experiments that 


the Briggs and Young correlation predicted the air heat transfer 
coefficients obtained from their experiementa! data for induced draft 
for six tube-row banks, but had to include correction factors given 
by Ward and Young (12) for tube banks having a number of rows 
less than six. 

Biery's method for predicting the Colburn j-modulus agreed 
with the Mirkovic and Briggs and Young correlations to within plus 
or minus 15 percent over the range of values of tube diameter, tube 
pitch, fin thickness, fin height, and fin spacing from which these 
correlations were based. However, these correlations deviated from 
Biery's method by as much as 40 percent at extreme values outside 
the range of data. All three correlations were also in agreement by 
predicting an increase in the Colburn j-~modulus with an increase in 
fin spacing or tube pitch or with a decrease in fin height or tube 
diameter. 

However, as fin thickness was increased, Biery's method and 
the Mirkovic correlation predicted an increase in the Colburn 
jJ-modulus, while the Briggs and Young correlation predicted a 
decrease. The experimental data taken by Mirkovic (7) presents an 
opportunity to examine variation of the Colburn j-modulus with fin 
thickness. In Mirkovic's run 12 and 13 all parameters were 
identical except for fin thickness. Mirkovic's Figure 2 shows for 


run 13 (the thicker fins), the Colburn j-modulus is higher than for 
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run 12. Therefore, an increase in the Colburn j-modulus' with 
increase in fin thickness seems to be the correct indication. 

In conclusion, all three correlations yield general agreement, 
except for variation of fin thickness. For this parameter two of the 
three correlations (Biery and Mirkovic) are in agreement both with 
each other and with a particular example from the Mirkovic data. 
Therefore, the choice of a correlation for ten the Colburn 
j-modulus of a given surface, rests between the Biery method and 
the Mirkovic correlation. In the final analysis, Biery's method was 
selected over the Mirkovic correlation. Biery's method is more 
general and based on a larger body of data encompassing a larger 
range of geometries. Further, it is not restricted to geometries in 
which the longitudinal distance between tube centers must be 
greater than one tube diameter as in the Mirkovic correlation. 
Therefore, the Biery method for predicting the Colburn j-modulus is 
used later for predicting the performance parameters of a given 


surface. 
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APPENGAIVS it 

The purpose of this Appendix is to compare Jameson's (5) cor- 
relation for pressure drop across triangular pitch finned tube banks 
with the data given by Kays and London (2). 

As noted in Chapter IV, to compare the friction factor used in 
the Jameson correlation, ar with the friction factor given in the 
data by Kays and London, a geometry factor must be included and 
the Reynolds number converted to a common diameter. With these 


changes the friction factor as defined by Kays and London (2) per 


equation (30) and as predicted by the Jameson correlation becomes: 


D., D, me eas (A23) 
fe eee OZ X_D nae 
ae. h 
Table All presents the percent deviation of the friction factor 
as predicted by the Jameson correlation, (f'), and the experimental 


value given by Kays and London (2) for the corresponding surface 


over the given range of Reynolds numbers, where 


Percent Deviation = fat X 100% (A24) 
Table All shows the friction factor predicted by the Jameson 


correlation (f') for surface 98D deviates from the experimental value 
by approximately 67%. Surface 98D is the one case of the surfaces 
listed in Table All which the minimum free flow area occurs in the 
diagonal spaces. Jameson (5) noted that his correlation did not 
accurately predict the friction factor for finned tubes in which the 


minimum free flow area occurs in the diagonal space. This 
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TABLE AI| 
Deviation of the Jameson (5) correlation for friction factor in 
predicting experimentally obtained value as presented in Kays and 


Wemdon (2). 


Surface*™ RE=1000 RE=4000 RE=7000 RE=10000 


92 ie On a7 +23. 87, +20 .0% +17.7% 
93 +17.3% +9.7%, +6 .7% +4 9%, 
94 +9 .9% ls cy SOno +2.4% 
96 +8 .9% +9 .8% +10.1% +10.47, 
97A Hal sows -0.1% -0.4% -0.5% 
97B -4,2% 231 184 =0.87, -0.2% 
98A S24 +0.7% +1.9% +2.7% 
98B -0.2% +0.3% +0.4% +0 .5%, 
98C +10.8% tl Oe Sy OO 4 +9 .9% 
98D +65.1% +67 .2% +68.1% +68 .6% 
98E 16.27 +7 .8% +8.4% +8 .9% 
99A +18.1% +17.7% +17.6% wll als 
99B +9.7%, +10.6% +11.0% +11.3% 


*Surface designation corresponds to the Kays and London (2) figure 


number. 





108 


conclusion appears to be supported by the results in Table All for 
surface 98D. Discounting this one surface, Table All shows the 
Jameson (5) correlation predicts the Kays and London (2) data to 
within approximately 20% of the mean value with a mean 7.8% above 
the Kays and London data and with a standard deviation of 7.3%. 
The Jameson correlation is used in Chapters V and VI to study the 
geometry effects of tube banks in which the minimum free flow 


area occurs in the spaces transverse to the air flow. 
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APPENDIX ey 

The purpose of this Appendix is to put together the equations 
and correlations introduced or derived in other Chapters 9 and 
Appendices in ane algorithm which may be programmed into a 
computer or calculator for predicting the performance parameters of 
a triangular pitch circular finned tube bank in crossflow given the 
geometry and the air and fin material properties. Hence, the only 
inputs required are tube diameter, tube pitch, fin thickness, fin 
height, fin spacing, and the air and fin material properties. 

The air and =“*fin material properties are entered as a 
"oroperty constant" et in the computation of fin effictency as 


outlined in Figure 5. 


Since: 
25 
_ 2h, (0) 
k § 
m 
mG: 
h N 2/3 USS) 
os 2/3 m2. Pr 
JO ° NST Pr @ ie 
cr 
where from (3a): 
RE p (A26) 
oe = D 
lt can be shown: 
2P_ j RE D 3 (A27) 
= c O . He) 
ml = Are ( f 
aue® 
where: 
aC, (A28) 
me - N 2/3, 
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FOr (Gmoemair at 90° F at atmospheric pressure ae Sel? 107¢ 
Tor aluminum fins and e —eOees on” hor COpber fims. 

The algorithm for predicting the performance parameters of a 
surface given the geometry and property constant is introduced in 
Figure A4&. This algorithm calculates the performance parameters for 
a triangular pitch smooth tube and finned tube with radial fins of 
rectangular profile only. An algorithm that would include tapered 
TINS, continuous fins, in-line tubes, and other various geometries 
could be written employing the methods outlined in this paper, but 
would require the applicable correlations and geometric relation- 
ships. It should be recalled that this algorithm is not accurate for 
predicting enhanced surface performance parameters for minimum 
Tree Tlow area occuring in the diagonal due to the inaccuracy of 
the Jameson (5) correlation for this case. However, the program 
should produce accurate results for both cases of minimum free flow 
Tor bare tubes. Therefore, the effect of varying tube pitch can be 
studied even for the case of minimum free flow occuring in the 
diagonal by using the performance parameters produced for smooth 
tubes. The program is written to produce performance parameters for 
bare tube Reynolds numbers (RE, |) of 10000, 70000, and 130000. 
Values for only two Reynolds numbers are required to plot the 
performance parameters for a surface since the plot is a straight 
line on a log-log scale. However, a third value insures plotting 
aeceurmacy. Ihe range for RE, of 10000 to 130000 was chosen since 
the typical value of RE, /RE is ¢ to 20, therefore, performance 


Darameters are produced for Reynolds mumber based on the 
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hydraulic diameter in the range of 500 to 18000. This is the 
approximate range most finned tube heat exchangers in crossflow 
operate and is the same range for which the Kays and London (2) 
data is given. 

The algorithm of Figure A&4 was keyed into a TI-59 calculator. 
This calculator program is listed in Figure A5. Due to the length of 
the algorithm, the calculator program is divided into two halves. 
The first half calculates the parameters for the equations in the 
algorithm up to label A (Steps 000 to 422), and stores these values 
in data registers for use by the second half of the program. After 
the first half of the program has been run, the second half of the 
program (Steps 000 to 632) can be read in by magnetic card and 
run to generate the performance parameters. Table AIlll lists the 
data registers in which the input and calculated parameters are 
stored. Where more than one parameter is listed indicates’ the 
corresponding data register is used for temporary storage of the 
indicated parameter while the program is running. 

Figure A6& shows a sample output generated by running the 
program of Figure A5. The program labels refer to the indicated 
parameter. The surface is defined in the output by a listing of the 
Tape parameters. It should be noted that the only input parameter 
not listed in the output its Poe In the parametric study of Chapter 
Vl the fin material was assumed to be copper and the fluid on the 
enhanced surface side was assumed to be dry air at Some for all 
cases. The property constant (es for these assumptions is 6.23 


ome: 
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The output of Figure A6 can be used to plot the performance 
parameters for the indicated surface. A comparison can then be 
made on the surface effectiveness as outlined in Chapter 2. Each 
surface generated in the parametric study of Chapter VI was 
compared on the basis of volume one side for the same pumping 
power and number of transfer units (Case C comparison of Chapter 


l1). Each surface was compared to a bare tube surface where X_ = 


- 
Xp = 2 and Dd, = 1". To standardize this comparison the 45 degree 
line (slope = 1 on log-log scale) was drawn from the point on the 


performance parameter plot of the reference bare tube surface where 


log (NTU/V = 2.0; 


| 

The 45 degree line described above and the plot of the 
performance parameters of a surface are straight lines. Therefore, 
it should be finally mentioned that as an alternative to plotting a 
performance parameter to make a volume comparison to another 
surface, linear equations may be written and solved simultaneously. 


This permits an analytical comparison which is faster if in the form 


of a program and also more accurate. 
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Figure A4: Algorithm for Prediction of Performance Parameters for 


Triangular Pitch Circular Finned Tube Surfaces and Bare 
Tube Surfaces. 
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